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0\ : Abstract 

We present a family of localized radiation modes in multilayered periodic media, where in-phase 

O . superposition of p-polarized waves leads to radiative confinement around the beam axis. Excitation 

• ^ . 

Qh' of surface plasmon polaritons yields an enhanced localization normally to the interfaces. On the 

o : 

c/3 ■ other hand, the spectral filtering induced by the presence of bandgaps allows to achieve transverse 

O 

c/3 ■ superresolution on the interfaces. Subwavelength beamwidths along an infinitely long distance 

[ might potentially be obtained. 
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I. INTRODUCTION 

Bessel beams [l-3| are a family of nondiffracting waves propagating in dielectric media 
along for instance the z-axis with a homogeneous spatial variation exp {i(3z) but transversally 
localized around its focus, say r = 0. For a zero-order Bessel beam, the transverse pattern 
is written in terms of the radially-symmetric function Jq {k±r). They are interpreted as 
a suitable superposition of plane waves all having a wavevector which projected onto the 
propagation axis yields the characteristic propagation constant /3 > of the Bessel beam 
4, |5| . A phase matching condition at focus is additionally satisfied in order to provide the 
best spatial localization of the wave in the vicinities of such an axis. The FWHM of the 
Bessel beam is 3.04 fcj^, which increases along with the propagation constant in accordance 
with the dispersion equation k'j_ + P'^ = n^k^, being n is the index of refraction of the medium 
and /cq = 27r/Ao the wavenumber in vacuum. Therefore the minimum beamwidth, roughly 
given by Ao/2?t,, is reached in the stationary- wave limit /3 = 0. 

This sort of solutions of the Helmholtz equation may be found also in stratified media if 
the unit vector (n) normal to the interfaces lies along the z-axis leading to normal incidence 
of the wavefield |6|-l9|. However, out-of-plane excitation cannot support the invariant prop- 
agation of Bessel beams 10|. For example, assuming the medium is periodic, the elements 
of any wave superposition are necessarily Bloch modes leading to a significant asymmetry 
in the resultant pattern. Provided the projection /3 of the pseudomoment along a direction 
perpendicular to the periodicity of the medium coincides for all Bloch components we may 
have a nondiffracting localized beam if, additionally, the phase matching condition is satis- 



fied 



With dielectric slabs, we encounter a higher localization of the field for decreasing 



values of the propagation constant (5. This is also reported in two-dimensional photonic 
crystals 12, 13 1, being observed not only under subdiffractive light propagation jl^ but 



much more generally. 

In this paper we study the formation of localized waves which propagation is resonantly 
sustained on a given interface from a ID periodic metallo-dielectric medium. The surface 
resonance arises in transverse-magnetic waves (H ± n) provided the sign of the dielectric 
constant changes abruptly at both sides of the interface. The excitation of such surface 
plasmons polaritons (SPPs) are attained at comparatively high values of /3, however leading 
to a subwavelength beam size. In spite of such an extreme wave localization the filamentation 
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of the excited field in tlie linear medium is ideally maintained for infinitely-long distances. 



II. IN-PLANE PROPAGATION OF DIFFRACTION-FREE BEAMS IN METALLO- 
DIELECTRIC LAYERED MEDIA 

Let us consider a monochromatic nondiffr acting beam propagating in a multilayered 
medium. The y-axis is set such that it is perpendicular to the surfaces separating the metallic 
media and the adjacent dielectric media, so that n = y. In Fig. [T]we show schematically the 
multilayered system. The width of a metallic slab is w, an element periodically replicated 
at a distance p along the y-axis. We also assume that beam propagation is directed along 
the 2;-axis, so that we may cast the electromagnetic fields as 

E(x, y, z, t) = e{x, y) exp {iPz — iujt) , (la) 
H(a:, ?/, 2;, t) = h{x, y) exp {iPz — iujt) , (lb) 

being u the frequency of the monochromatic radiation. The homogeneity of the wave field 
in the coordinate z is explicitly parametrized in terms of the propagation constant (3. More 
specifically, we study TM waves where H exists only onto planes parallel to the metal- 
dielectric interfaces, and therefore which component hy vanishes. Later on we additionally 
consider wave confinement around the origin {x,y) = (0,0) in any given transverse plane. 

The Maxwell's equations provide some relations between the transverse fields e and h. 
The electric field e may be derived from h by means of the equation V x H = —iueoeE, 
where e{y) is the relative dielectric constant of the foliar structure. Thus we have 

= i (weoe)"^ dyh^, (2a) 
Cy = i {ueoey^ {if3h^ - d^h^) , (2b) 
62 = -i {ueoeY^ dyh^. (2c) 

From the equation V (/xoH) = we also find 

= i(3~'^d^h^, (3) 

so that h;j. is the scalar wavefield from which we may describe the nondiffr acting beam 
unambiguously. Let us point out that alternate routes for the description of electromagnetic 



diffraction-free beams may be found elsewhere 
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The two-dimensional Helmholtz equation 



(92 + 9^ + 0; Woe-/?') K = Q. 



(4) 



may be of help in order to find any localized solution of the field h^. However it is preferable 



an alternative approach. Based on the Bloch theorem [18 



we may describe the propagating 



wavefield as a superposition of Bloch modes having the form 

/oo 
aKhxiy) exp {it^^x + iKy) dk^ 
-00 



(5) 



where is the (real) component of a wave vector along the x-axis, and K is the so-called 
Block wave number, which generally is multivalued for a given kx- Moreover, hxiy + p) = 
hxiy) is a periodic function, which is normalized as hxiO) = 1 for convenience, and ax is 
simply a field amplitude. 

Assuming that sources are sufficiently far from the focal region of the localized beam, the 
role of evanescent waves [l^ is negligible in our analysis and, therefore, they are disregarded 
setting ax = if Im(i^') 7^ 0. In this sense we also neglect material losses imposing that the 
relative dielectric constant of the medium is a real parameter of the problem [Im(e) = 0] . We 
point out that this is not strictly true in a realistic problem; plasma-like media with e < are 
necessarily dispersive since energy density considerations lead to the condition d{ijje)/doj > 0, 
and the presence of dispersion in general signifies dissipation of energy [20]. However, it is 
possible to neglect the absorption in a transparency window where the imaginary part of e 
is very small in comparison with its real part. 

In the case of the stratified medium shown in Fig. [1], the periodic function is conveniently 
written as the addition of two independent terms, hx = + hj^. Using the column vector 



h 



K 



'■K 



(6) 



we may write 21] 



h.K{y) = exp {-iKy) Ta{y - ya)^a, for y e Rc 



(7) 



where a is an integer that refers to a unique slab. The domains = {ya-i,ya] where 
yo = 0, yi = w and ya+2 = Va + P', for instance, Ra is associated with a region where the 
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medium is metallic when a is an odd number. In Eq. ([7]), the translation matrix 



T.(l/) 



exp {iky^y) 







(8) 







exp {-ikyaU) 



and the column vector Vq, = [a^, ha\ giving some constant amplitudes. Note that is 
independent of given in Eq. (|5]) and they should not be confused. Also, 



where ei = and e2 = are the relative dielectric constants of the metallic medium and 
the dielectric material, respectively. For other slabs we use the recurrence relation ea+2 = ^a] 
here ec at a = does not refers to its value in vacuum but and it is not necessarily the 
unity. In Eq. the real-valued wave vector onto a plane parallel to the metallo-dielectric 
interface has a squared modulus fcy = fc^ + Finally since material losses are neglected, 
the dielectric constant of the plasma-like material < and that of the insulator > 
are also real-valued parameters. 

We also take advantage of spatial symmetries from the layered system in Fig. [H In 
particular, we use that the relative dielectric constant is a real and even function e{—y') = 
e{y'), being y' = y — ya-1/2 a spatial coordinate centered in the middle point of the domain 
Ra, ya-1/2 = iVa + ya-i) /2. Givcu a wavcfield hx{x,y') that satisfies Eq. (jl]), the fields 
hl{x,y') and hx{x, —y') are also solutions of the wave equation. From the series expansion 
of Eq. ([5]), therefore, we infer that hx^—y') and h*j^{y') are periodic functions both associated 
with the same Bloch (real) number —K, so that there should be a constant of proportionality 
between them. Provided ip is the argument of hxiy' = 0) we find that 



As a consequence hx is hermitian under the normalization hxiy' = 0) = 1. However this 
is not necessarily true when hxiv = 0) = 1. In the later case it yields for instance that the 
amplitude of the periodic function evaluated on every surface y = ya reaches the unit. 

The continuity at the interfaces of the wavefields h;j. given in Eq. i^, h^, and [that 
is continuity of e'^dyh^ as inferred from Eqs. (l2a|) and fl2cl) ]. together with the periodicity 
of hx determine the waveform of each Bloch mode and its dispersion behaviour. Using the 




(9) 



hK{-y')e^v{-i^) = [hK{y')Qw{-i'^)] 



(10) 
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transfer matrix elements for layered media 2l| we may obtain in a rather simple way the 
dispersion equation. Let us give the corresponding analysis briefly. Experienced readers 
may ignore the rest of the present section thus going to Sec. [1111 directly. 



A. Dispersion equation and isofrequency curves 



Eq. ([7]) provides the complex waveform of h^; appropriately. However, in order to sim- 
plify the application of continuity conditions of the fields, we also establish an equivalent 
definition, 

hK{y) = exp i-iKy) T^{y - ya-iWa (H) 
at y E Ra, where the new column vector = [a^, b'^]. Note that 

= PaV„, (12) 

the matrix = Ta{ya-i — ya) representing a translation between the boundaries of 
Thus the continuity of the p-polarized fields at ?/ = ?/„ holds when 



where the transmission matrix 



■\jyoL ""ya 



(13) 



(14) 



Consequently, the field amplitudes corresponding to the slab labelled as (a + 1) may be 
determined by means of those amplitudes of the adjacent slab a as v^+i = Pq^^iD^^^Dq,Vq,. 
On other hand, imposing the periodicity condition /ix(ya+2) = h^iya) we also find 



v„+2 = exp (iKp) Vo 



(15) 



given as a direct application of the Bloch theorem. As a consequence, considering the 
continuity of the fields at y = ya given in ( fT3l) and that at y = ya+i, Eq- ( |T5l) leads to the 
matrix equation 

M^Vq = exp (iKp) v„, (16) 

where Eq. (|T2l) is also used at (a + 1) and (a + 2). The computation of the matrix Mq, leads 

to 



Ca 



(17) 
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Periodicity in the medium renders Pq,+2 = Pa and = Da which yield further sim- 

phfications in the determination of Mq,. We leave the reader the long but straightforward 
evaluation of the elements of the matrix Mq,. In general M„_(_i 7^ Mq,; however M„+2 = M^j 
and, more importantly, the matrix Mq is unimodular [det(MQ,) = 1] and its trace is inde- 
pendent of a. 

In agreement with Eq. (fT6!) . Vq, is an eigenvector of the matrix Mq, providing the eigen- 
value exp {iKp) . Consequently the determinant of the matrix Mq — exp {iKp) I vanishes, 
where I is the identity matrix. The resultant equation may be cast simply as cos (Kp) = 
{Aa + Da) /2, where the unimodularity property of the matrix Mq, AaDa — BaCa = 1, 
should be employed. Thus we have derived the dispersion equation, which explicitly reads 



cos {Kp) = cos[kyd{p — w)] cos{kyrnw) (18) 

sm[kyd{p - w)] sm{kymw). 



2,kyi^kyfYi6(j^6lJl 

As expected, the dispersion equation (ITS]) is independent of the integer a. Our interest is 
focused in real values of K, for which the right-hand side of Eq. ( 1181) should be comprised 
between —1 and 1. 

In Fig. [2] we plot the real- valued K solutions of equation ( !T8|) for a multilayer medium 
composed of thin metallic films of width w = 50 nm, separated a distance p = 450 nm, and 
embedded in a dielectric medium of = 2.25. At a frequency to = 3.427 fs"^ (wavelength 
Aq = 550 nm in the vacuum) our plasma-like medium has a relative dielectric constant 
em = —15.0. A large bandgap in the interval k\\p G (0,4.571) is followed by some two others 
in (5.126, 7.899) and (8.301, 8.503), together with the evanescent-wave regime at k\\p > 8.737. 
As a consequence, real values of (3 are limited by /3max = 8.737/p in our example {f3, 



19.42 fim~^). This boundary value is reached at K = ivr/p rather than at a zero value 22 |. 
Let us point out that, if the superposition shown in Eq. is such that k\\ is higher than 
the cut-off frequency 

/3c = V^dio/c, (19) 

occurring if k\\p > 7.711, the wavefields are all of evanescent nature in the metal and in 
the dielectric; here Pc = 17.14 fim~^. The cut-off frequency is plotted with a horizontal 
red line in Fig. O Interestingly, these evanescent waves may arise at values of k\\p such 
that Eq. (ITS!) gives a Bloch wave number with Im(i^) = 0. Thus, the evanescent waves 
are resonantly coupled leading to propagating Bloch-type constituents of the diffraction-free 



wavefield hx{x,y). Nondiffracting beams with a propagation constant (3 G {Pc, Pmax) are 
wave fields of tliis kind. 

In Fig. [3] we map different contours of isofrequency /3 in tlie k^K plane based on the 
graphical representation of the dispersion equation given in Fig. [21 Thus we consider real 
values of K again. The first band is represented in Fig. [21(a) providing a surface of maxima at 
kx = Q and K = ±tt/p reaching /3max- Simple closed curves around these points are found for 
decreasing values of the propagation constant provided (3ci < (3 < f3max, where (3cip = 8.503 
coincides with the upper boundary of the first bandgap. At lower f3, spatial frequencies 
around k^ = cannot excite propagating Bloch modes leading to open (isofrequency) curves. 
This gap is removed in the second sheet of the dispersion equation, which arises for /3 < /3c2 
and it is shown in Fig. [H^b); obviously f3c2P = 8.301 corresponds to the upper limit of the 
second band of allowed spatial frequencies /cy. In this case, the closed curves are centered 
at the origin of the plane k^K (and at multiples of 271 /p along the i^-axis). Similarly to 
Fig. [31(a), the isofrequency curves are closed when /3 takes values belonging to the (second) 
band of allowed frequencies, here /3c3 < < (3c2 where (3c3P = 7.899. Out of this interval, 
the isofrequency curves are open again. 



B. Bloch components of the wavefield 

According to Eq. (fTSl) . once the amplitudes vq and vi are given for all Bloch wavenumbers 
K, the periodic function hx and thus the wavefield are estimated unambiguously. In fact 
we have shown that vi = P^^D^^DoVq, therefore the degrees of freedom are reduced to 
the amplitudes vq. With the normalization /^/^(O) = 1 we additionally impose that the 
summation of the elements of vq is left to unity. As a result, the amplitudes associated with 
the Bloch modes in the dielectric slab at Rq = {—p + w,0] yield 

1 

exp (iKp) + Bq- Aq 

where Aq and Bq are matrix elements of Mq. 
We point out that the denominator in Eq. (1201 

kym ( C 

kyd 



Bo 

exp (iKp) — Aq 



(20) 



) vanishes if 



(21) 
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which in principle might lead to a singular behavior of vq. In this case we also have Bq = 
0. Eq. fl2T]) represents the condition for the existence of a surface resonance in a metal- 
dielectric interface 23|. Therefore kyd would be purely imaginary, which is associated with 



the existence of evanescent waves in the dielectric films. In addition Mq would be a diagonal 
matrix with 

±K= ^ky^ -{}-^) ^y'^ (22) 
providing the eigenvector vq as either the unit vector [1,0] or [0, 1]. Moreover, Im(i^') ^ 
except when K vanishes, that occurs if additionally 

= (23) 

ed w 

leading to Mq = I. The latter particular case might bring new insights onto our basic 
analysis however it deserves a thorough study out of here. 

The periodic function hx is depicted in Fig. HJfor the multilayered setup with parameters 
given in the caption of Fig. |2]and at different values of k\\ . In all cases, the complex amplitude 
{IikI is an even function whereas its argument is an odd function respect to the center of 
the metallic film, as shown in Eq. (fTOl) . provided it is set off to zero. Peaks are formed 
at metal-dielectric interfaces and valleys at the center of the slabs; this rule holds in the 
purely-evanescent regime when k\\ > Pc, however the maximum of intensity is reached at 
the center of the dielectrics if fcy belongs to the third band [see Fig. 111(e)]. On other hand 
we observe that phases of the wavefunction have a characteristic variation at different 
spectral bands. For instance, phases are ~ vr rad at the right side of the central metallic 
slab y = w in the top band; however, this phase turns to ~ in the second band. We show 
below that this is of relevance in the formation of an on-axis focus. 



III. FOCUS GENERATION 

Let us establish some favourable conditions for the formation of a focus along the z-axis. 
Using {x,y) = (0,0) in Eq. we obtain the wavefield amplitude = / (^Rdkx at 
the origin as a summation of the amplitudes ax corresponding to different Bloch modes. 
This may be interpreted as an interference of Bloch-type individuals. If the phase of their 
amplitudes are manipulated in order to have the same value leading to in-phase waves, the 
oscillatory superposition yields the highest intensity achievable. Excepting a few particular 
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conditions, it cannot be found a point other than the origin from the xy plane where such 
a phase matching holds. As a consequence, a strong localization of the nondiffr acting beam 
is expected to occur around the 2;-axis, such a line unquestionably constituting a focus. 

At this point of our analysis it is interesting to review the concept of focus wavefields 
in a system such as the uniform dielectric medium, where nondiffracting solutions of the 
wave equation are well known. Obviously we are speaking of Bessel beams ij. Such a 
system would result from setting = in the layered medium of Fig. [1] to remove the 
metallic films from the dielectric host. In this case Eq. ( ITSjl leads trivially to the solutions 
K± = ±kyd. Also hxiy) = 1 in the whole xy plane. The phase matching condition at the 
origin is also observed over a radially-symmetric transverse pattern if in addition ax is in 
direct proportion to \d(f)/dk^\ = \ — 1/K\, where is the polar angle in the kr^K plane. Thus 
inserting 

if \k^\ < ^/Pl 



V/3|-(fei+/32)' 

0, otherwise 



(3' 



(24) 



into Eq. ([5]) yields the Bessel wavefield hx = 27!-Jo{k±r) provided /3 < f3c [see Eq. (fT9|) ]. being 



k^ = .M^ 



(25) 



the transverse wave number and r the radial spatial coordinate. See Ref. 10| for further 



lm{K) = 0, aK 







details. We point out that excitations of the type ax = [i 
otherwis e] a lso provide paraxial Bessel beams in periodic media and in anisotropic 
crystals [24 1. 



A field spectrum like that of Eq. ( I24l) may be experimentally attained using an opaque 



screen, with a centered extremely-thin transparent annulus, placed at the front foca 



mirrors 
elements 



a perfect lens j4|. Highly-efficient approaches may be found usin g co nical lenses 



25 



plane of 



)ry-Perot interferometers 



26| and 



leaky screens 29l-l3l| . and diffractive optical 



32|-|34||. Using such devices as external sources in our system would excite the 



required diffraction-free wavefields in the layered medium. For simplicity we assume that 
the spectral strength of such nondiffracting beams has a form following Eq. ([21 

1 



for \kx\ < 7, 



(26) 



and dK = if \kx\ > 7- 
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Let us consider that /3 and 7 may be tuned at convenience. Thus (3 <k\\< ^/W 

+ 7^ so 

that we might apply any spectral band (i.e. ordinates in Fig. |2]) arbitrarily. In our model it 



is reasonable to think of \fW^^r^ denoting the cut-off frequency /3c = 27rA/ed/Ao associated 
with the uniform external medium of dielectric constant q where the source field is driven. 
In a great number of examples given below /3c = f^max yielding q = 2.89; in this medium 
the wavelength = Ao/y^ = 324 nm and the minimum spot size of the zero-order Bessel 
beam is (Aa;)^^^ = 2.253/7 = US nm where the maximum value of 7 = /3c at /3 = 0. In 
other hand, our model neglects filtering (apodizing) effects and aberrations on the wave 



fields induced at the boundaries of the system [10|. However, conventional techniques for 
its compensation might be employed in order to find a good agreement with our results. 
Otherwise the theory remains valid leaving the appropriate estimation of the spectrum ax- 



IV. SUB WAVELENGTH TRANSVERSE PATTERNS: NUMERICAL RESULTS 

To illustrate the focus generation along the z-axis, we perform a numerical simulation in 
the periodic layered medium of p = 450 nm (and = 2.25) with metallic films of w = 50 nm 
and era = —15.0. In Fig. [5] we show the field intensity l/i^P for a nondiffracting beam of 
propagation constant /3 = 18.90 yum^^ (/3 = Pci), which spatial spectrum is given by Eq. ( 126|) 
with 7 = 4.462 fim~^ {'jp = 2.008). Thus the maximum value of fcy reaches a value coinciding 
with f3max, i-e. /3ci < k\\ < f3max corresponding to the first band of allowed frequencies shown 
in Fig. El The isofrequency curve K = K{kx) is also depicted in the inset of the figure, 
where excited frequencies kx are shaded in blue. The intensity is maximum at the origin 
but the pattern exhibits no radial symmetry. The field distribution along the abscissa, 

hx{x.,y = Q) = 2 I cos (7a; cos 0) (i0 = 27r Jo(7x), (27) 

resembles that of the source field that would propagate in the homogeneous dielectric medium 
of Q = 2.89. We point out that some discrepancies arise in Eq. (127]) if some spatial fre- 
quencies kx are excited but, because of the presence of a gap, they do not contribute to the 
resulting propagating wavefield; this case will be treated ahead. From Eq. ( 1271) we derive 
that the FWHM of the intensity peak along this direction is inversely proportional to 7 
following Aa, = 2.253/7. Ill our case A^^ = 505 nm, which is above the wavelength in the 
dielectric host A^ = Aq/i/q (= 367 nm). However, the behavior along the ordinate is signif- 
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icantly different. The most attractive feature of the wavefield is its high locahzation in the 
metal-dielectric interfaces, leading to fast decays when moving away from the surfaces and 
thus forming wedge-like shapes. This is in agreement with the patterns shown in Figs. SJ^a) 
and (b) contributing to the integral at x = 0. Although the highest peak is attained at 
y = 0, a large one also arises on the other side of the central metallic film, y = w. There, the 
Bloch modes hxiy) exp (iKy) are as strong as in focus and they interfere nearly in-phase 
(dephase < 0.297r rad) giving a secondary focus. Ignoring this sidelobe, the FWHM of the 
figure is Aj, = 45.99 nm, well below A^. 

Apart from surface resonances, radial asymmetry from focus is also attributed to a high 
concentration of light along certain directions in the xy plane. This directional enhancement 



of the radiated power may be explained by means of an effect coined as photon focusing 



35|, 



a term which first was used onto the strong anisotropy of heat flux in crystalline solids and 



that takes other names in optics like self-collimation 36|], self-guiding [STj, and sub diffr active 



propagation 



The estimation of the value(s) of the azimuthal coordinate 6 where photon focusing is 
manifested is based on the stationary-phase principle. Firstly note that the phase factor 
exp {ik^x + iKy) of Eq. (|5]) varies rapidly when moving away from the origin, r — t- oo, being 
{x,y) = r {cos 6, sin 6). This phase term also introduces strong oscillations for running k^ 
except in the vicinities of those spatial frequencies k^g satisfying x + K'^y = 0, where K' = 
dK/ dkx and the subindex s stands for the value given at k^s- Setting K'^ = tan (ps we have the 
solutions (ps = 9±7i/2; consequently the tangent of the curve K{kx) at the stationary points 
kxs is normal to the vector position of the observation point (x, y). Moreover, the asymptotic 
behavior of Eq. ([5]) depends exclusively on those stationary points k^s] substituting K by 
Ks in every term of the integrand excepting the phase factor for which we use a quadratic 
expansion K ^ Ks + K'^ik^ — k^g) + K'^ik^ — k^sY l'^, finally yields 

Kix,y)^ ^ axhKXy) exp {ik^sX + iK^y) iRsiv), (28) 



where 



/oo 
exp [iK'^ih, - kxsfy/2\ dk, 
-oo 




i-sign{K':y) . (29) 
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In this simple analysis we have also assumed that ax has a well behavior and it does not 
present discontinuities. In the far field zone l/i^P decreases inversely proportional to the 
radial coordinate r, excepting those directions 6 where stationary points satisfy K'^ = 0. 
In these cases the dispersion curve is fiat leading to a significantly slow attenuation of the 
radiated power, at least much slower than r"^. 

In Fig. ini^a) we plot the solutions k^s of the equation K'^ = at different values of p. 
Considering that the dispersion equation f lTSj) may be written as cos{Kp) = /(fc^), the 
solutions of the equation K'J = are also the frequencies k^g satisfying 

fs (1 - fs) + f'sfs = 0. (30) 

As shown in Fig. |3]^a), the isofrequency curves for f3ci < f3 < (3max mimic ellipses that lack 
of fiat sections so that K" = has no real solutions. At /3 = /3ci, the isofrequency curve 
represented also in the inset of Fig. |5] is fiat at the origin of the plane k^K taking a shape 
in X. In this case K.,) = (0, 0) for which K'^ = ±0.604 and 0, = ±31 deg. See Fig. ^b) 
showing the estimated (positive) values of (pg also for other propagation constants, /3. Our 
analysis conclude that the self-guiding phenomenon is expected at angles 6 = ±59 deg (and 
6 = ±121 deg), in agreement with the transverse pattern shown in Fig. |5l When Pc3 < P < 
Pel, solving Eq. (15CT]) leads to increasing values of \kxs\ and therefore increasing values of 
as shown in Figs. EJ^a) and (b) respectively. Interestingly, kxs = becomes again a solution 
if /3 = Pc3 (here giving = ±n/p) for which (pg = ±42 deg; also k^g = ±3.32/p gives 
Kg = where (pg = ±73 deg. Moreover, a minimum of four values of k^g may be found for 
P < Pc3 in which the zero-diffraction condition K" = holds. We may conclude the general 
rule that a larger number of solutions comes out when (positive) P decreases reaching new 
bandgaps for the frequency k\\ as shown in Fig. |5J 

Let us analyze the intensity patterns of diffraction-free beams under the presence of 
bandgaps [Im(ii') 7^ 0] within the source spectral window |A;^| < 7. In Fig. [7]we plot Ih^l"^ of 
wavefields with propagation constant P = Pcs corresponding to the lower limit of the second 
band. In Fig. [Tl^a) ■yp = 2.552 that yields {k\\)^^^ = Pc2- This means that the source excites 
entirely the second band of allowed spatial frequencies with no gaps; therefore Eq. ( !27l) is 
valid and A^; = 2.253/7 (= 397 nm). Along the y-axis, the narrow peak around the origin 
has a FWHM Ay = 62.52 nm and again is accompanied by a high sidelobe at the other 
side of the metallic film. The formation of this secondary focus is again the constructive 
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interference (nearly in-phase) of the Bloch modes from the second band. On other hand, 
self-colhmation is attributed exclusively to the stationary points {kxs, Kg) = (0, ±7r/p) giving 
9 = ±48 deg {(ps = ±42 deg as seen above). 

Setting 7p = 3.148 as shown in Fig. Wijo), the bandgap (3c2 < k\\ < (3^ is encountered, 
that is frequencies satisfying 7ci < \kx\ < are frustratedly excited (being 7cip = 2.552 
and 7c2P = 3.148), so that the optical system behaves like a low-pass filter. In this case 

hx{x, y = 0) = 27r Jo(7x) - h^ciiix), (31) 

where 

Kcl hx)=2 COS (7X COS 0) = ^ J2m M (32) 

being < 0ci < vr/2 such that 7cos0ci = Id, 

Ci-^=4(-l)- ""^';;^'^"\ form^0, (33) 

and Cq^^^ = 40ci- We point out that (0) = for m 7^ so that in the vicinities 
of the focal point we have ^ {2tt — A(j)ci) Joil^) ■ However the central peak stretches 
slightly, specifically = 459 nm, whereas sidelobes are altered significantly as expected. 
Moreover, FWHM in the direction of the periodicity Ay = 61.58 nm and self-colhmation 
angles 9 = ±48 deg are also maintained is spite of bandgaps. 

Finally, arriving at 7p = 3.734 we are able to excite every allowed spatial frequency in the 
first band since {k\\)^^^ = (3max- The intensity Ikx]"^ is depicted in Fig. [7](c) showing notable 
differences from the cases previously analyzed. For instance, the FWHM along the x-axis is 
A^. = 277 nm, which is close to the value we would obtain if we ignore the bandgap. In this 
case Ax is lower than the radiation wavelength in the dielectric A^^ = 367 nm; however this 
subwavelength size is still higher than Ad/2 in the same order than a regular Bessel beam. 
Here we may evaluate the wavefield by 

hx{x, y = 0) = 27rJo(7x) - Kciiix) + hxc2{ix), (34) 

where 7cos0c2 = 7c2- Sidelobes are also strongly attenuated in all directions; even at those 
associated with self-collimation. Here 9 = ±48 deg (0^ = ±42 deg) and also 6' = ±17 deg 
{(j)s = ±73 deg); as a consequence, sharing the power radiated in this large number of 
directions leads to the weakening of the photon-focusing effect. More importantly, the high 
sidelobe appearing previously in the ?/-axis seems to be wiped out completely. This effect 
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may be explained considering that Bloch components from the first band and those from 
the second band are roughly out of phase ai y = w [see Figs. |l](b) and (d)] so that they 
interfere destructively in Eq. Additionally, the central peak is unaltered in practical 
terms, giving A^^ = 54 nm. 

V. HYBRID CONSTRUCTION OF DIFFRACTION-FREE WAVES 

In the numerical simulations given above we have shown that beamsizes along the x-axis 
are larger than the diffraction limit Xd/^ attained by quasi- stationary Bessel beams propa- 
gating in the medium of dielectric constant e^, whereas Ay is clearly subwavelength. Control 
over the wave pattern and thus over its FWHM in the x direction is exercised by the spec- 
trum of spatial frequencies k^: the higher bandwidth the lower A^;. From k'^ < P^^x ~ P'^ 
estimating the spectral domain of the wavefield we conclude that decreasing /3 leads to a 
widening of the spectrum. However, if /3 < /3c the plane-wave components in the dielec- 
tric material propagate homogeneously so that they are not necessarily coupled resonantly 
around the metal-dielectric surfaces. In principle, this effect might modify significantly the 
localization of the field in the direction of the periodicity. Let us clarify these aspects. 

The contour plot shown in Fig. |8] corresponds to a diffraction-free beam of propagation 
constant /3 = 10.157 fim~^ [Pp = 4.571) coinciding with the lower limit of k\\ for the third 
band (see Fig. [2]). The transverse frequency •yp = 7.446 guarantees that all the three bands 
of allowed frequencies are excited; this is a high value since if it were generated in a uniform 
medium with q = 2.89 would give a Bessel beam of width A^, = 136 nm. Fig. |H] shows that 
the intensity is distributed mainly around the interface y = 0. A narrow peak on the 
focus is produced exhibiting a width A^^ = 132 nm, which is slightly lower than that just 
given above. Also strong sidelobes arise in the vicinity of focus. This suggests that bandgaps 
inherent in a photonic crystal provide a mechanism to achieve a superresolving effect going 
beyond the diffract ion- induced Bessel limit. The origin of this superresolving effect is not 
the excitation of SPPs but an enhancement of high spatial frequencies j39 |. 

On other hand, the field distribution in the y-axis demonstrates a subwavelength focus 
of FWHM Ay = 44 nm. Sidelobes on the interfaces are accompanied with others peaks in 
the middle of the dielectric slabs. This is not surprising since Bloch components of the third 
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band shown in Fig. |H^e) and contributing in the expansion (|5]) have such a behavior. This 
effect is in fact more pronounced along directions oriented with polar angles 6* = 0^ ± 7r/2 
being the stationary point 0^ = =1=41 deg {0 = ±49 deg) associated with the self-coUimation 
effect. 

Although superresolution along the x-axis is modest in comparison with that obtained 
in the y-axis, we wonder if beamwidths might surpass the diffraction limit in all directions. 
As resolved in Sec. IIIIl for our examples where Pmax = 19.42 fim~^ and Aq = 550 nm this 
diffraction limit reads A^in = 116 nm. In Fig. [9]we represent the intensity distribution 
of a nondiffracting beam excited with a transverse spatial frequency 7 = 17.20 yum^^ giving 
a propagation constant (3 = 9.00 /im"^. In this case the central peak is anamorphic and 
characterized by FWHMs A^. = 109 nm and A^ = 46 nm all below the diffraction limit. The 
bandgap- induced filtering also leads to intense sidelobes along the interface y = 0, however 
they decay as going beyond the focus. Moreover, the peak intensity of these secondary foci 
cannot exceed that intensity of the main focus. 

VI. CONCLUDING REMARKS 

We have demonstrated that diffraction-free beams propagating in structured media com- 
posed of alternating layers of positive and negative e may reach beamsizes surpassing the 
diffraction limit. For an excitation of the form given in ( |26ll . we recall that this limit is 
attained when the transverse spatial frequency is maximum, jmax = y/^^/c, leading to 
stationary (/3 = 0) Bessel waves with a central spot of FWHM Amin = 2.253/7maa'- In 
the layered system, two different mechanisms lead to the superresolving effect. A bandpass 
filtering due to the existence of gaps in the spatial spectrum of modifies the response 
of the system transversally providing a narrow peak along the x-axis with moderate gain 
(A^ < Amin) and high sidelobes (secondary foci). In the direction of the periodicity, however, 
superresolution is carried out by the formation of surface resonances in the metal-dielectric 
interfaces leading to fast decays out of these planes. In our numerical simulations. Ay can be 
as low as a third of A^m- Additionally, dephase of Bloch constituents belonging to different 
bands contributes to the growth control of secondary foci in nearby surfaces. 

In principle, subwavelength beamwidths along an infinitely long distance might poten- 
tially be obtained. However practical limitations in the geometry of the system leads to the 
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generation of invariant wavefields along with a finite axial distance. The estimation of such 
imits is similar for example to that corresponding to Bessel beams propagating in free space 
2I, and therefore they are not discussed here. However we point out the decisive role of 
the structured medium in order to achieve the subwavelength beamsizes: these wavefields 
cannot be sustained in free space for ranges longer than a wavelength. Relevant aspects in 



near-field subwave^ 
and co-workers 



ength nondiffracting beams are exhaustively investigated by Kukhlevsky 
4o|-|43|. 

Finally, material losses are neglected in our analysis. Since Im(e) 7^ leads to wave 
damping in an unbounded dissipative system, absorption in the medium should be accounted 
for assuming a limited number of slabs. We have performed numerical simulations showing 
that the superresolving effect along the y-axis is practically maintained for a small absorption 
coefficient of the metallic medium and a moderate number of layers, whereas the peak slightly 
gets wider along the x-axis. Even assuming that the focus is placed at an interface near the 
center of the system, boundary effects and also tunnelling effects cannot be ignored. For 
instance, setting = — 15 + i0.15 we may excite a diffraction-free beam with a propagation 
constant (3 = 9 fim~^ and a central peak of A^^ = 139 nm and Ay = 43 nm in a medium 
composed of 20 unit cells. The deterioration in resolution power along the x-axis is mainly 
caused by a selective increment in Im(i^) within the allowed spectral bands, an effect that 
is specially dramatic in the highest band. We have estimated that the central peak is also 
enlarged to a lesser extent by increasing the number of layers. Under these unfavourable 
conditions the beamsize is still surpassing the diffraction limit along the direction of the 
periodicity. On other hand, a significant attenuation of the focal intensity is also evident. 
In our numerical simulation, in-focus intensity decreases down to 18% (down to 6%) from 
its original value if Im (e^) = 0.15 and the number of unit cells is 10 (20). A manuscript 
containing an extended analysis with further details is in preparation. 



17 



ACKNOWLEDGMENTS 



This research was funded by Ministerio de Ciencia e Innovacion (MICIIN) under the 
project TEC2009-11635. 



[1] J. Durnin, J. Opt. Soc. Am. A 4, 651 (1987). 

[2] J. Durnin, J. J. Miceli, and J. H. Eberly, Phys. Rev. Lett. 58, 1499 (1987). 

[3] Localized Waves, edited by H. E. Hernandez-Figueroa, M. Zamboni-Rached, and E. Recami 

(John Wiley & Sons, 2008). 
[4] G. Indebctouw, J. Opt. Soc. Am. A 6, 150 (1989). 

[5] J. Fagerholm, A. T. Friberg, J. Huttunen, D. P. Morgan, and M. M. Salomaa, Phys. Rev. E 
54, 4347 (1996). 

[6] D. Mugnai, A. Ranfagni, and R. Ruggeri, Phys. Rev. Lett. 84, 4830 (2000). 

[7] S. Longhi, K. Janner, and P. Laporta, J. Opt. B 6, 477 (2004). 

[8] A. V. Novitsky and D. V. Novitsky J. Phys. A: Math. Gen. 39, 5227 (2006). 

[9] C. J. Zapata-Rodn'guez, M. A. Porras, and J. J. Miret, J. Opt. Soc. Am. A 25, 2758 (2008). 
[10] C. J. Zapata-Rodrfguez and J. J. Miret, arXiv:0910.5642vl [physics.optics](2009). 
[11] J. J. Miret and C. J. Zapata-Rodn'guez, J. Opt. Soc. Am. B 25, 1 (2008). 
[12] D. N. Christodouhdes, F. Lederer, and Y. Silberberg, Nature 424, 817 (2003). 
[13] O. Manela, M. Segev, and D. N. Christodouhdes, Opt. Lett. 30, 2611 (2005). 
[14] K. Stahunas and R. Herrero, Phys. Rev. E 73, 016601 (2006). 
[15] Z. Bouchal and M. Ohvik, J. Mod. Opt. 45, 1555 (1995). 
[16] Z. Bouchal, J. Bajer, and M. Bertolotti, J. Opt. Soc. Am. A 15, 2172 (1998). 
[17] P. Paakkoncn, J. Tcrvo, P. Vahimaa, J. Turuncn, and F. Gori, Opt. Express 10, 949 (2002). 
[18] J. D. Joannopoulos, S. G. Johnson, J. N. Winn, and R. D. Meade, Photonic crystals. Molding 

the flow of light (Princeton University Press, 2008). 
[19] S. A. Ramakrishna and O. J. F. Martin, Opt. Lett. 30, 2626 (2004). 

[20] L. D. Landau, E. M. Lifshitz, and L. P. Pitaevskii, Electrodynamics of continuous media 

(Butterworth-Heinenann, Oxford, England, 1984). 
[21] P. Yeh, Optical Waves in Layered Media (Wiley, New York, 1988). 



18 



[22] V. Kuzmiak and A. A. Maradudin, Phys. Rev. B 55, 7427 (1997). 

[23] H. Raether, Surface plasmons on smooth and rough surfaces and on gratings (Springer- Verlag, 
Berlin, 1988). 

[24] A. Ciattoni and C. Palma, Opt. Commun. 224, 175 (2003). 

[25] R. Grunwald, U. Griebner, F. Tschirschwitz, E. T. J. Nibbering, T. Elsaesser, V. Kebbel, 

H.-J. Hartmann, and W. Jiiptner, Opt. Lett. 25, 981 (2000). 
[26] C. J. Zapata-Rodrfguez and A. Sanchez-Losa, J. Opt. Soc. Am. A 23, 3016 (2006). 
[27] K. B. Kuntz, B. Braverman, S. H. Youn, M. Lobino, E. M. Pessina, and A. I. Lvovsky, Phys. 

Rev. A 79, 043802 (2009). 
[28] Z. L. Horvath, M. Erdelyi, G. Szabo, Z. Bor, F. K. Tittel, and J. R. Cavallaro, J. Opt. Soc. 

Am. A 14, 3009 (1997). 
[29] S. Holm, IEEE Trans. Ultrason., Ferroelect., Freq. Contr. 45, 712 (1998). 
[30] K. Reivelt and P. Saari, Phys. Rev. E 65, 046622 (2002). 
[31] K. Reivelt and P. Saari, Phys. Rev. E 66, 056611 (2002). 
[32] A. Vasara, J. Turunen, and A. T. Friberg, J. Opt. Soc. Am. A 6, 1748 (1989). 
[33] J. Amako, D. Sawaki, and E. Fujii, J. Opt. Soc. Am. B 20, 2562 (2003). 
[34] Z. Li, K. B. Alici, H. Caglayan, and E. Ozbay, Phys. Rev. Lett. 102, 143901 (2009). 
[35] D. N. Chigrin, Phys. Rev. E 70, 056611 (2004). 

[36] H. Kosaka, T. Kawashima, A. Tomita, M. Notomi, T. Tamamura, T. Sato, and S. Kawakami, 

Appl. Phys. Lett. 74, 12121214 (1999). 
[37] D. N. Chigrin, S. Enoch, C. M. Sotomayor-Torres, and G. Tayeb, Opt. Express 11, 1203 

(2003) . 

[38] K. Staliunas, C. Scrrat, R. Herrero, C. Cojocaru, and J. Trull, Phys. Rev. E 74, 016605 (2006). 
[39] M. Martfnez-Corral, P. Andres, C. J. Zapata-Rodrfguez, and M. Kowalczyk, Opt. Commun. 
165, 267 (1999). 

[40] S. V. Kukhlevsky, M. Mechler, L. Csapo, and K. Janssens, Phys. Lett. A 319, 439 (2003). 
[41] S. V. Kukhlevsky, M. Mechler, L. Csapo, K. Janssens, and O. Samek, Phys. Rev. B 70, 195428 

(2004) . 

[42] S. V. Kukhlevsky and M. Mechler, Opt. Commun. 231, 35 (2004). 

[43] G. Nyitray, V. Mathew, and S. V. Kukhlevsky, Opt. Commun. 281, 1082 (2008). 



19 



LIST OF FIGURES 



1 Schematic geometry of the planar-nanolayer-based medium |23 

2 Dispersion equation at a; = 3.4 x 10^^ rad/s for a periodic media as represented 
in Fig. [T]with w = 5 x 10"^ m and p = 4.5 x 10^^ m. The red hues marks the 
boundary of homogeneous- wave and evanescent- wave regimes in the dielectric, 

fc||P = 7.711 I24 

3 Isofrequency curves at different propagation constants f3c < (3 < Pmax are 
shown in (a) for the first sheet and in (b) for the second sheet of the dispersion 
curve. Contour hues are labelled following the normalization /3p |25 

4 Behavior of the complex amplitude (left column) and argument (right column) 
of the wavefunction hx at different values of k\\ and positive K corresponding 

to the first band (top), second band (middle) and third band (bottom) |26 

5 Contour plot of the intensity \hx\'^ in the xy plane corresponding to a local- 
ized diffraction-free beam of normalized propagation constant /3p = 8.503 and 
transverse frequency 7p = 2.008 leading to = 8.737. Intensity distri- 
butions along the coordinate axes are shown at the top and left sides. Inset: 
Isofrequency curve where shaded region corresponds to the excited spatial 
bandwidth |27 

6 (a) Solutions k^. of the equation (130|) . K'J = 0, associated with different prop- 
agation constant /3 of the wavefields, in the layered structure of Fig. [H (b) 
Azimuthal angle in the plane k^K that corresponds to each solution of the 
aforementioned self-guiding condition |28 

7 Transverse intensity \hx\'^{x,y) for diffraction-free beams of Pp = 7.899 and 
spatial frequency (a) ■yp = 2.552, (b) •yp = 3.148, and (c) yp = 3.734. Insets: 
Isofrequency curve at f3p = 7.899. The blue shade indicates the spectral win- 
dow in kx under excitation. Light blue regions refer to gap-induced frustrated 
excitation 29 



20 



Intensity \hx\'^ in the transverse xy plane for a nondiffr acting beam of /3 = 
10.16 yum~^ and 7 = 16.55 /im^^. Excitation of Bloch modes with \kx\ < 
5.15 /im~^ leads to nonevanescent wavelets in the dielectric slabs, so that 
the resulting hybrid wavefield from Eq. ^ combines Bloch constituents of 

different nature 

Transverse pattern of a nondiffracting beam with /3 = 9.00 /im~^. The 
beamwidth of the main central focus has A^; = 109 nm and Ay = 46 nm, 
going beyond the diffraction limit A^m = 116 nm 



21 



FIGURES 



22 



FIG. 1. Schematic geometry of the planar-nanolayer-based medium. 
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FIG. 2. Dispersion equation at uj = 3.4 x 10^^ rad/s for a periodic media as represented in Fig. [1] 
with w = 5x 10^^ m and p = 4.5 x 10^^ m. The red hnes marks the boundary of homogeneous- wave 
and evanescent-wave regimes in the dielectric, fcyp = 7.711. 
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(a) (b) 




FIG. 3. Isofrequency curves at different propagation constants I3c < f3 < Pmax are shown in (a) for 
the first sheet and in (b) for the second sheet of the dispersion curve. Contour lines are labelled 
following the normalization jip 
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FIG. 4. Behavior of the complex ampUtude (left column) and argument (right column) of the 
wavefunction hx at different values of fcy and positive K corresponding to the first band (top), 
second band (middle) and third band (bottom). 
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FIG. 5. Contour plot of the intensity \hx\'^ in the xy plane corresponding to a localized diffraction- 
free beam of normalized propagation constant (3p = 8.503 and transverse frequency 7p = 2.008 
leading to = 8.737. Intensity distributions along the coordinate axes are shown at the 

top and left sides. Inset: Isofrequency curve where shaded region corresponds to the excited spatial 
bandwidth. 
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(a) (b) 




FIG. 6. (a) Solutions kx of the equation ()30p . K'J = 0, associated with different propagation 
constant /3 of the wavefields, in the layered structure of Fig. [TJ (b) Azimuthal angle 4> in the plane 
kxK that corresponds to each solution of the aforementioned self-guiding condition. 
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FIG. 7. Transverse intensity \hx\'^{x,y) for diffraction-free beams of /3p = 7.899 and spatial fre- 
quency (a) 7p = 2.552, (b) jp = 3.148, and (c) 7p = 3.734. Insets: Isofrequency curve at 
/3p = 7.899. The blue shade indicates the spectral window in kx under excitation. Light blue 
regions refer to gap-induced frustrated excitation. 
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FIG. 8. Intensity \h x\ in the transverse xy plane for a nondifFracting beam of /3 — 10.16 /ini 
and 7 = 16.55 fim^^. Excitation of Bloch modes with \kx\ < 5.15 fim^^ leads to nonevanescent 
wavelets in the dielectric slabs, so that the resulting hybrid wavefield from Eq. ([5]) combines Bloch 
constituents of different nature. 
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FIG. 9. Transverse pattern of a nondiffracting beam with (3 = 9.00 //m . The beamwidth of 
the main central focus has A^^ = 109 nm and Aj^ = 46 nm, going beyond the diffraction hmit 
Amm = 116 nm. 
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